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Abstract 

A high statistics sample of photoproduced charm particles from the FOCUS (E831) 
experiment at Fermilab has been used to search for CP violation in the Cabibbo 
suppressed decay modes D + — > K~ K + ir + , D° — > K~K + and D° — > 7r~7r + . 
We have measured the following CP asymmetry parameters: Acp(K~ K + ir + ) = 
+0.006 ± 0.011 ± 0.005, A C p{K-K+) = -0.001 ± 0.022 ± 0.015 and A C p(tt'tt + ) = 
+0.048 ± 0.039 ± 0.025 where the first error is statistical and the second error is 
systematic. These asymmetries are consistent with zero with smaller errors than 
previous measurements. 
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1 Introduction 



CP violation occurs if the decay rate for a particle differs from the decay rate 
of its CP-conjugate particle [1]. CP violation, which in the Standard Model 
(SM) is a consequence of a complex amplitude in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, has been observed to date only in the neutral-kaon 
system. In charm meson decays (as well as in K and B decays), two classes 
of CP violation exist: indirect and direct. In neutral-charm- meson decays, in- 
direct CP violation may arise due to D° — D mixing. In the case of direct 
violation, CP violating effects occur in a decay process only if the decay am- 
plitude is the sum of two different parts, whose phases are made of a weak 
(CKM) and a strong contribution due to final state interactions (FSI) [2]. The 
weak contributions to the phases change sign when going to the CP-conjugate 
process, while the strong ones do not. In singly Cabibbo-suppressed D decays, 
penguin terms in the effective Hamiltonian may provide the different phases 
of the two weak amplitudes. 

Compared to the strange and bottom sectors, the SM predictions of CP vi- 
olation for charm decays are much smaller [2-5], making the charm sector a 
good place to test the SM and to look for evidence of new physics. In the SM, 
direct CP violating asymmetries in D decays are predicted to be largest in 
singly Cabibbo-suppressed decays, at most 10~ 3 , and non-existent in Cabibbo- 
favored and doubly Cabibbo-suppressed decays [1]. However, a CP asymme- 
try could occur in the decay modes D — > K s mi due to interference between 
Cabibbo-favored and doubly Cabibbo-suppressed decays. 
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Experimentally, one looks at the Cabibbo-suppressed decay modes which have 
the largest combination of branching fraction and detection efficiency For this 
reason we select the all-charged decay modes D + — > K~K + n + , D° — > K~K + , 
and D° — > it~it + (throughout this paper the charge conjugate state is implied, 
unless otherwise noted). 

In D decays the charged D is self-tagging and the neutral D is tagged as either 
a D° or a D by using the sign of the bachelor pion in the D*^ decay. 

Before searching for a CP asymmetry we must account for differences, at the 
production level, between D and D in photoproduction (the hadronization 
process, in the presence of remnant quarks from the nucleon, gives rise to 
production asymmetries [6] ). This is done by normalizing to the Cabibbo- 
favored modes D° — > K~n + and D + — > K~-n + -n + , with the additional benefit 
that most of the corrections due to inefficiencies cancel out, reducing system- 
atic uncertainties. An implicit assumption is that there is no measurable CP 
violation in the Cabibbo-favored decays. 

The CP asymmetry can be written as: 

V(D) - rj(D) 



i-CP 



rj(D)+rj(D) 



(1) 



where r] is (considering for example the decay mode D° — > K K + ): 

N(D° -»• K-K+) e(D° -> R-tt + ) 



T)(D) 



N(D° -> K-TI+) e(D° -> K-K+) 



where N(D° — > K K + ) is the number of reconstructed candidate decays and 
e(D° — > K~K + ) is the efficiency obtained from Monte Carlo simulations. 

The CP asymmetry parameter measures the direct CP asymmetry in the case 
of D + and the combined direct and indirect CP asymmetries for D° [7] . 

The name FOCUS stands for Photoproduction of Charm with an Upgraded 
Spectrometer with a lexical license. The word "upgrade" refers to the upgrade 
of the E687 (the predecessor experiment) spectrometer [8]. 

Charmed particles were produced by the interaction of high energy photons, 
obtained by means of bremsstrahlung of electron and positron beams (with 
typically 300 GeV endpoint energy), with a beryllium oxide target. The mean 
energy of the photon beam was approximately 180 GeV. The data were col- 
lected at Fermilab during the 1996-97 fixed-target run. More than 6.3 x 10 9 
triggers were collected from which more than 1 million charmed particles have 
been reconstructed. 
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The particles from the interaction are detected in a large-aperture magnetic 
spectrometer with excellent vertex measurement, particle identification and 
calorimetric capabilities. The vertex detector consists of two systems of sil- 
icon microvertex detectors. The upstream system consists of 4 planes inter- 
leaved with the experimental target, while the downstream system consists 
of 12 planes of microstrips arranged in three views. These detectors provide 
high resolution separation of primary (production) and secondary (decay) ver- 
tices with an average proper time resolution of 30 fs for 2-track vertices. The 
momentum of the charged particles is determined by measuring their deflec- 
tions in two analysis magnets of opposite polarity with five stations of multi- 
wire proportional chambers. Kaons and pions in the D-meson final states are 
well separated up to 60 GeV/c of momentum using three multicell threshold 
Cerenkov counters. 



2 Selection criteria 

The final states are selected using a candidate driven vertex algorithm [8]. A 
secondary vertex is formed from the reconstructed tracks and the momentum 
vector of the D candidate is used as a seed to intersect the other tracks in 
the event to find the primary vertex. Once the production and decay ver- 
tices are determined, the distance i between them and the relative error <jg 
are computed. Cuts on the Ijan ratio are applied to extract the D signals 
from the prompt background. The topological configuration of the event is 
evaluated with four tests: the primary and secondary vertex confidence levels 
(minimum values of 1% were required) and two measures of vertex isolation, 
a no point-back isolation and a secondary vertex isolation. The no point-back 
isolation cut requires that the maximum confidence level for a candidate-/} 
daughter track to form a vertex with the tracks from the primary vertex be 
less than a certain threshold. The secondary vertex isolation cut requires that 
the maximum confidence level for another track to form a vertex with the D 
candidate be less than a certain threshold. The analyses differ mainly in the 
way the particle identification is handled and, less importantly, in the way 
the vertex cuts are applied. To minimize the systematic error we use identical 
vertex cuts on the signal and normalizing modes. 

In the D + — > K~K + tt + analysis, we require £/<Jg > 10, the no point-back 
isolation must be less than 20%, and the secondary vertex isolation less than 
0.1%. The vertices (primary and secondary) have to lie inside a fiducial vol- 
wm^J the D momentum must be in the range 25 < P < 250 GeV/c (a 

16 The reason for this cut lies in the presence of a trigger counter just upstream of 
the second microstrip device, therefore we defined the fiducial volume as the target 
region between the first slab of the experimental target and this trigger counter. 
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Fig. 1. (a) K~K +/ k + invariant mass distribution, (b) K + K~ir~ invariant mass 
distribution, (c) K~tt + tt + invariant mass distribution, (d) K + tt~tt~ invariant mass 
distribution. The fits (solid curves) are described in the text and the numbers quoted 
are the yields. 

very loose cut) and the primary vertex must be formed with at least two 
reconstructed tracks, in addition to the seed track. The Cerenkov particle 
identification cuts used in FOCUS are based on likelihood ratios between the 
various stable particle identification hypotheses. These likelihoods are com- 
puted for a given track from the observed firing response (on or off) of all 
cells within the track's {(3 = 1) Cerenkov cone for each of our three Cerenkov 
counters. The product of all firing probabilities for all cells within the three 
Cerenkov cones produces a x 2 -like variable Wi = —2 In (Likelihood) where i 
ranges over the electron, pion, kaon and proton hypotheses (see reference [9] 
for more details). We require W n — Wk > f for the tracks reconstructed as 
a kaon, and a pion consistency cut for the pion tracks. The pion consistency 
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cut requires that no particle hypothesis is favored over the pion hypothesis 
with a AW = W n — W m i n exceeding 2. To remove contamination from the 
— > K~K + n + signal due to Cerenkov misidentified background from the 
decay mode D + — > K~it + -n + , we employ an anti-reflection cut which rejects 
candidates which, when reconstructed as K~it + -n + , are consistent with the 
D + hypothesis (we also reject events whose K~K + mass exceeds 1.84 GeV/c 2 
in order to exclude background due to D* + — > D°tt + — > (K~ K + )tt + ). This 
cut has no effect in the vicinity of the D + — > K~K + tt + signal peak. 

In Fig.l, the invariant mass plots obtained with this set of cuts for the decay 
modes D + — > K~K + tt + , D~ — > K + K~tt~ , and the normalizing decays D + — > 
K~7t + tt + and D~ — > fT + 7r~7r~ are shown. In the D + — > A'~7r + 7r + analysis 
there is an additional cut on the D° mass formed by a kaon and a pion to 
remove the D* + — > D°7r + — > (X ~7r + )7r + decay chain. The KKn invariant 
mass distributions are fit with a Gaussian for the D + signal, a second Gaussian 
for the Df signal, and a quadratic polynomial for the background. A binned 
maximum likelihood fit finds 6860±110 D + -> K-R+1T+ and 7355±112 /J- -> 
K + K~ir~ events. The fit for the normalizing modes (fit with a Gaussian plus 
a linear polynomial) gives 68607±282 D+ -> fs:-7r + 7r + and 73710 ±292 D~ -> 
A' + 7r _ 7r _ events. 

In the D° — > K~K + analysis, the sign of the bachelor pion in the D*^ decay 
chain — > D°(D is used to identify the neutral Z) as either a £>° 

or a D°. We require that the mass difference between the D° and the D* mass 
be within 4 MeV/c 2 with respect to the nominal mass difference [10]. We use 
i/o-f > 8, while the no point-back isolation and the secondary vertex isolation 
cuts are unnecessary because the D* tag sufficiently reduces the background. 
All the other cuts, except the anti-reflection cuts, are the same as those used 
in the D + — > K~K +r n + analysis. 

In Fig. 2, we show the invariant mass plots, obtained with this set of cuts, for 
the decay modes D° — > K~K + , D — > K + K~, and the normalizing decays 
D° — > K~n + and D — >• fT + 7r~. The peak in the KK invariant mass plots 
at f=s 1.95 GeV/c 2 is due to the reflection of the D° — * K~tt + mode when 
the pion is misidentified as a kaon. A Monte Carlo simulation of this reflection 
reproduces the shape observed in the data. Consequently, the K~ K + invariant 
mass distributions are fit with a Gaussian for the D° signal, a function obtained 
by smoothing the reflection peak (only the shape of this reflection peak is 
modeled by our Monte Carlo simulation, the amplitude of this peak is given 
by a free parameter of the fit) and a quadratic polynomial for the background. 
From a binned maximum likelihood fit we find 1623 ± 47 D° — >■ K~K + and 
1707 ± 53 D — > K + K~ events. The fit for the normalizing modes (fit with 
a Gaussian plus a linear polynomial and excluding the low mass region to 
avoid possible contamination due to other charm hadronic decays involving 
an additional vr° ) gives 18501±144 D° -> and 19633±149 if -> K+n- 
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Fig. 2. (a) K~K + invariant mass distribution, (b) K + K~ invariant mass distribu- 
tion, (c) K~ir + invariant mass distribution, (d) K + n~ invariant mass distribution. 
The fits (solid curves) are described in the text and the numbers quoted are the 
yields. 



events. 



The D° — > 7r~7r + analysis is identical to the D° — > K~K + analysis with the 
exception of the Cerenkov identification. In order to reduce the large reflection 
peak from the D° — > K~n + mode, a tight Cerenkov identification requirement 
(Wk — W w > 1) for the pion tracks is implemented. 

In Fig. 3 the invariant mass plots for the decay modes D° — > 7r~7r + and D° — > 
7r + 7r _ are shown. As in the D° — > K~K + case, the peak at ~ 1.75GeV/c 2 
is due to the reflection of the D° — > K~tt + mode. Again our Monte Carlo 
simulation reproduces the shape observed in the data. The 7v~7i + invariant 
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Fig. 3. (a) ir~ir + invariant mass distribution, (b) tt + tt~ invariant mass distribution. 
The fits (solid curves) are described in the text and the numbers quoted are the 
yields. 

mass distributions are fit with a Gaussian for the D° signal, a function ob- 
tained by smoothing the reflection peak, and a quadratic polynomial for the 
background. A binned maximum likelihood fit gives 606 ±31 D° — > 7r~7r + and 
571 ± 33 D — > 7r + 7r~ events. 



3 CP asymmetry measurements 

Because the efficiency is strongly dependent on the D momentum, it is nec- 
essary to verify that the observed momentum spectrum is reproduced by the 
Monte Carlo simulation. A mismatch could generate a false asymmetry. Fig. 4 
shows the D° momentum for the decay mode D° — > K~n + for real data (points 
with errors) and Monte Carlo data (histogram). The D° momentum spectrum 
is obtained (both in real data and Monte Carlo data) by subtracting the events 
contained in the sideband regions from the events in the signal region. The 
two histograms are normalized by area. The agreement is good. 

The measured asymmetries are reported in table 1 along with a comparison to 
previous measurements. The analysis of the decay mode D + — > K~K + n + is 
complicated by the possibility of intermediate resonant states, such as K*°K + 
and 07r + . We think that in this case a Dalitz plot analysis is the appropriate 
tool to investigate CP asymmmetry effects. This topic will be the subject of 
a future paper. 

Our asymmetry measurements have been tested by modifying each of the 
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Fig. 4. D momentum for the decay mode D — > i^T 7r + . The points with errors 
correspond to the real data, the histogram to the Monte Carlo data. 

vertex cuts individually; the results are always consistent. 
Table 1 



Measured CP asymmetries (xlO 2 ). 



Experiment 


D + -> K-K+7T+ D° -> K~K+ D° -► 7^+ 


E687 [11] 
CLEO II [12] 
E791 [13,14] 
This measurement 


-3.1 ±6.8 +2.4 ±8.4 
±8.0 ±6.1 

-1.4 ±2.9 -1.0 ±4.9 ±1.2 -4.9 ±7.8 ±3.0 
±0.6 ±1.1 ±0.5 -0.1 ±2.2 ±1.5 ±4.8 ±3.9 ±2.5 



The systematic errors on these measurements reflect uncertainties in recon- 
struction efficiency, Cerenkov particle identification, and hadronic absorption 
of secondaries in the target and spectrometer materials. The estimates were 
obtained by splitting our data into independent samples depending on D mo- 
mentum and the different periods in which the data were collected. The main 
reason for this time depedence is the insertion of the upstream silicon system 
(in the target region) during the 1997 fixed-target run period. A technique 
modeled after the S-factor method from the Particle Data Group [10] was 
used to separate true systematic variations from statistical fluctuations. The 
asymmetry is evaluated for each of the statistically independent subsamples 
and a scaled variance is calculated; the split sample variance is defined as the 
difference between the reported statistical variance and the scaled variance if 
the scaled variance exceeds the statistical variance. For all decays studied, the 
scaled variance is less than the statistical variance, so it does not contribute 
to the systematic uncertainty. The evaluation of systematic effects related to 
different fit procedures is performed on the whole sample. The asymmetry is 
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Table 2 



Contributions to the systematic uncertainty for the measured asymmetries. 



Source 


A CP (K-K + -K + ) 


A CP {K-K+) 


A C p(ir 7T+) 


Split sample 


no contribution 


no contribution 


no contribution 


Fit variant 


0.0016 


0.0144 


0.0243 


Monte Carlo statistics 


0.0050 


0.0054 


0.0055 


Total systematic error 


0.0053 


0.0154 


0.0249 



calculated using various fit conditions (these include the choice of the estima- 
tor, the background shape and, in the case of the D° analysis, the shape of 
the reflection peak) and the sample variance is used because the fit variants 
are all a priori likely. To obtain the final systematic error, the variance from 
the different fitting procedures and a further contribution, due to the uncer- 
tainties in the efficiency calculation and finite statistics of the Monte Carlo 
events, are then added in quadrature. Table 2 shows the contribution of each 
of these sources to the total systematic uncertainty for the CP asymmetry 
measurements. 

The measured asymmetries are consistent with zero within the errors. 



4 Summary and Conclusions 

We have searched for CP violation in the Cabibbo-suppressed decay modes 
D + — > K~K + n + , D° — > K~K + and D° — > it~it + using a high statistics sam- 
ple of photoproduced charm particles from the FOCUS (E831) experiment 
at Fermilab. We have measured the following CP asymmetry parameters: 
A CP (K-K + TT + ) = +0.006 ± 0.011 ± 0.005, A CP (K~ K + ) = -0.001 ± 0.022 ± 
0.015 and A CP {w + ) = +0.048 ± 0.039 ± 0.025 where the first error is sta- 
tistical and the second error is systematic. 

These asymmetries are consistent with zero and represent a substantial im- 
provement over previous measurements. 
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